%% ===========================������˵����=============================
%ǰ���ļ�------
%��������------
%�������------
%�������------
%% ===========================���������塿=============================
function [d0,Ac,L0,k]=MCiniRandom(VarMat,coeOfVar,svm,...
    usvm,index,simMCNum)
sigmaMat=VarMat.*coeOfVar;
varNum=length(svm); 
mnMat=m_n((1:simMCNum)',varNum);  
Md0=normrnd(VarMat(1),sigmaMat(1),[simMCNum,1]);
MAc=normrnd(VarMat(2),sigmaMat(2),[simMCNum,1]);
ML0=normrnd(VarMat(3),sigmaMat(3),[simMCNum,1]);
Mk=normrnd(VarMat(4),sigmaMat(4),[simMCNum,1]);
Mcell={Md0,MAc,ML0,Mk};
for ak47=1:varNum  
    Mcell{svm(ak47)}=Mcell{svm(ak47)}(mnMat(:,ak47));   
end
for m4a1=1:(4-varNum)  
    Mcell{usvm(m4a1)}=Mcell{usvm(m4a1)}(mnMat(:,1));
end
Md0=Mcell{1};
MAc=Mcell{2};
ML0=Mcell{3};
Mk=Mcell{4};
d0=Md0(index);
Ac=MAc(index);
L0=ML0(index);
k=Mk(index);
